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Branching fraction measurements using B-meson decays to Kg-K ■k~ are presented. These mea- 
surements were obtained by analyzing a data sample of 88.9 million T(4S) BB decays collected 
with the BABMi detector at the SLAC PEP-II asymmetric-energy B factory. Using a maximum 
likelihood fit, the following branching fraction results were obtained: B{B^ K^-k'^'k~) — (43.7 ± 
3.8 ± 3.4) X 10"®, B{B° K*+-k-) = (12.9 ± 2.4 ± 1.4) x 10"'^ and 6(5" KI-k-)-k+ = 

(42.7 ± 2.1 ± 2.2) X 10"''. The CP violating charge asymmetry Ak'-k for the decay B° K'+tt' 
was measured to be Ak'-^ = 0.23±0.18l:l^;!^^. For all these measurements the first error is statistical 
and the second is systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Three-body decays of the B meson tend to be domi- 
nated by intermediate quasi two-body charmed particles 
with the charmless resonant and non-resonant contribu- 
tion being small. Nevertheless, these charmless decays 
prove to be important in furthering our understanding 
of the weak interaction and complex quark couplings de- 
scribed by the Cabibbo-Kobayashi-Maskawa matrix ele- 
ments JJ. 

The i3-meson decay to K'^'!t~^tt~ can proceed via many 
interesting charmless resonances which we can probe for 
CP violation, such as foK° 0, p°K'^ and K*+Tr~ . A 
limit on the sum of their branching fractions can be 
obtained by measuring the inclusive charmless branch- 
ing fraction of B'^ — s- K'^tt^tt". This measurement has 
been performed previously by the CLEO 3] and Belle JJ| 
experiments. For the mode B'^ K*~^'!t~ the branch- 
ing fraction can be measured directly with the available 
BABAR data sample. 

Branching fraction and asymmetry measurements of 
charmless B decays can also be used to test the accuracy 
of QCD factorization models 5]. In particular there are 
factorization models that predict CP asymmetries in the 
decay B° K*+tt~ ^. The decay 5° ^ K*+tt- is self 
tagged (the charge of the kaon reflects the flavor of the 
_B-meson), so the CP asymmetry can be defined as: 



K' + Tv- 



(1) 



In this paper the branching fractions of B'-^ 



K°n+7T-, B° 



K 



and B° D-{-^ K°tt- 



are presented, where charge conjugate decays are also 
implied. The procedure used selection criteria requiring 
events with a reconstructed K'^TT^Tr~ final state. In the 
case of B{B'^ /^^tt+tt"), the total charmless contri- 
bution to the Dalitz plot was measured (with charmed 
and charmonium resonances removed), including con- 
tributions from resonant charmless sub-structure. For 
the decays B° K*+tt- and ^ D-{-^ K°Tr-)7:+, 
the analysis was restricted to the region of the K'^tt^tt^ 
Dahtz plot consistent with K*~^{-^ -f^^s""^) ^"^^ 
i^sTT") decays respectively. Finally, the Ak*tt value for 
the decay B'^ K*~^tt~ , which was first measured by 
CLEO .7] was extracted. 

The data used in this analysis were collected at 
the PEP-II asymmetric-energy e~^e~ storage ring with 



the BABAR detector |g. The BABAR detector consists 
of a double-sided five-layer silicon tracker, a 40-layer 
drift chamber, a Cherenkov detector, an electromagnetic 
calorimeter and a magnet with instrumented flux re- 
turn. The data sample has an integrated luminosity of 
81.9 fb^^ collected at the T(4S) resonance, which corre- 
sponds to (88.9 ± 1.0) X 10^ BB pairs. It was assumed 
that the T(4S) decays equally to neutral and charged B- 
meson pairs. In addition, 9.6 fb~^ of data collected at 
40 MeV below the T(4S) resonance were used for back- 
ground studies. 

Candidate _B-mesons were reconstructed from two 
tracks and a K'^, where the K*^ was reconstructed from 
7r+7r~ candidates. Each of the two tracks that were not 
generated by the K'^ were required to have at least 12 
hits in the drift chamber, a transverse momentum greater 
than 100 MeV/c and to be consistent with originating 
from the beam-spot. These tracks were selected as pious 
using energy loss {dE / dx) measured in the tracking sys- 
tem, the number of photons measured by the Cherenkov 
detector and their corresponding Cherenkov angle. Fur- 
thermore, the tracks were also required to fail the elec- 
tron selection based on dE/dx information, the ratio of 
energy in the calorimeter to momentum in the drift cham- 
ber and the shape of the signal in the calorimeter. The 
prerequisites imposed on K'^ candidates were for the re- 
constructed mass to be within 15 MeV/c^ of the nominal 
K'^ mass , a decay vertex separated from the B'^ decay 
vertex by at least five standard deviations and a cosine 
of the angle between the line joining the B and K'^ decay 
vertices and the K*^ momentum to be greater than 0.999. 

To characterize signal events, two kinematic and one 
event shape variable were used. The first kinematic vari- 
able AE, is the difference between the ccnter-of-mass 
(CM) energy of the B-candidate and y/s/2, where ^/s is 
the total CM energy. The second is the beam-energy- 
substituted mass Toes = \/ {s/2 + Pj • /Ef — p^, 
where ps is the B momentum and (Bi,pi) is the four- 
momentum of the T(4S) in the laboratory frame. Us- 
ing these two kinematic variables, candidates had to 
be in the range |AB| < 0.1 GeV and 5.22 < toes < 
5.29GeV/c^. The event shape variable is a Fisher dis- 
criminant {T) ^3 ■ The T variable was constructed from 
a linear combination of the cosine of the angle between 
the B-candidate momentum and the beam axis, the co- 
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sine of the angle between the B-candidate thrust axis 
and the beam axis, and the energy flow of the rest of the 
event into each of nine contiguous, concentric, 10° cones 
around the thrust axis of the reconstructed B 

Continuum quark production (e+e~ — > qq where q — 
u,d,s,c) was by far the dominant source of background. 
This was suppressed using another event-shape variable 
which was the cosine of the angle Ot between the thrust 
axis of the selected S-candidate and the thrust axis of 
the rest of the event. For continuum background, the 
distribution of | cos 9t \ is strongly peaked towards unity 
whereas the distribution is flat for signal events. There- 
fore, the relative amount of continuum background was 
reduced by requiring that all candidates fulfill the crite- 
rion I cos6't| < 0.9. 

Simulated Monte Carlo (MC) events were used to 
study background from other B-meson decays. The 
largest potential S-background was seen to come from 
quasi two-body decays including charmonium mesons 
such as J/ipK'^, XcqK^ and ■i/;(2S')ifg where the char- 
monium meson decays to n'^fJ-^ which are misidenti- 
fied as pions or where they decay directly to tt+tt". 
These background events were removed by vetoing re- 
constructed TT+TT" masses consistent with 3.04 < r7i^+^- 
< 3.17 GeV/c2, 3.32 < m^+^- < 3.53 GeV/c^ and 
3.60 < mTT+TT- < 3.78 GeV/c^, identifying the J/'tp, XcO 
and ^{2S) mesons respectively. Additionally, in order 
to measure the charmless branching fraction of the de- 
cay B° K°Tr+Tr-, B° D' {-^ K^sTt')tt+ events 
were removed by vetoing events with a reconstructed 
ii'gTr invariant mass consistent with 1.83 < rrij^o^ < 
1.90 GeV/c^. Monte Carlo simulation showed that 21 
± 3 i?'^ ^ D^{-^ i^37r^)7r+ background events still re- 
mained. These events had a reconstructed D~ mass out- 
side the veto as a result of using the wrong K^^ or 7r+ 
which was incorrectly selected from the other B decay 
in the event. When selecting B^ K*^7r^ or B° 
D~{-^ K^'K~)'K^ candidates, the additional cuts 0.79 < 
m^o^ < 0.99 GeV/c^ and 1.85 < m^o^ < 1.89 GeV/c^ 
were applied respectively to the reconstructed m^o^ in- 
variant mass. After the above selection criteria were ap- 
plied, a small proportion of events for all decays under 
study had more than one candidate which satisfied the se- 
lection criteria. For these events, one candidate alone was 
selected by choosing the candidate whose cos^r value 
was closest to 0. In a signal MC study, this selects the 
true signal candidate more than 75% of the time. 

After all cuts, the largest remaining _B-background 
to B° K°TT+Tr- was the 4-body decay B° r)'K° 
with r]' p° (770)7 and p° TT+TT which contributes 
22 ± 6 events. For the B" K*+ti- and B° 
D^{-^ Xs7r")7r+ channels, the background contribution 
was small and came from modes that can interfere by 
decaying to a K'^tt^tt^ final state such as foK'^ and 
p^K^. In addition the K*^'k~ and _D+7r~ modes are 
backgrounds to each other. Furthermore, there was the 



non-resonant if^Tr+Tr" background contribution to the 
resonant signal. Along with selection efhciencies obtained 
from MC, using available information on exclusive mea- 
surements or by fitting to regions in the Dalitz plot, 
upper limits or branching fractions for these modes were 
obtained to estimate their background contributions. 

In order to extract the signal event yield for the channel 
under study, an unbinned extended maximum likelihood 
fit was used. The likehhood function for TV candidates is: 

CM \ ^ / *^ \ 

where i, j and I are integers, M is the number of hypothe- 
ses (signal, continuum background and i?-background) , 
Pi (a, Xj) is a probability density function (PDF) with 
the parameters a depending on three variables (x) toes, 
AE, and J^, and ni is the number of events for each hy- 
pothesis determined by maximizing the likelihood func- 
tion. The PDF is a product Pi{d, Xj) = Pi{am-ES^ wes) • 
Pi{ctAE, AE) ■ Pi{ajr^T). Correlations between these 
variables were small for signal and continuum back- 
ground hypotheses. However for S-background, cor- 
relations were observed between toes and AE, which 
were taken into account by forming a 2-dimensional 
PDF for these variables. The parameters of the sig- 
nal and B-background PDFs were determined from MC. 
The continuum background parameters were allowed to 
vary in the fit, to help reduce systematic effects from 
this dominant event type. Upper sideband data de- 
fined to be in the region 0.1 < AE < 0.3 GeV and 
5.22 < Toes < 5.29 GeV/c^ was used to model the contin- 
uum background PDFs. For the toes PDFs a Gaussian 
distribution was used for signal and a threshold func- 
tion |0| was used for continuum. For the AE PDFs a 
sum of two Gaussian distributions with the same means 
was used for the signal and a first order polynomial was 
used for the continuum background. Finally, for the 
PDFs, a sum of two Gaussian distributions with distinct 
means and widths was used for signal and an asymmetric 
Gaussian which has different widths above and below the 
modal value was used for continuum background. In the 
case of i?-background parameterizations, signal-like or 
continuum-like PDFs were used depending on the char- 
acteristics of the background. With more than 400 signal 
events and typically a one-to-one signal to background 
ratio in the total number of B^ D^{-^ K^Itt^)^^ 
candidates, it was possible also to vary the signal PDF 
parameters in the fit for this mode. This enabled un- 
certainties and corrections due to MC, to be calculated 
and apphed to the B° K^n+n- and B° K*+tt- 
analyses. Figure ^ shows the fitted projections of the 
maximum likelihood fit to B'^ D^{-^ K'^tt^)tt^ can- 
didates in ruEs, and J- containing 472 ± 24 signal 
and 455 ± 23 background candidates. 
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ture above 1 GeV/c^ which is the region where higher 
kaon resonances can occur. 



5.29 
niEs (GeV/c') 



2.0 3.0 
Fisher 

FIG. 1: Maximum likelihood fit projections of mss (a), Ai5 
(b) and T (c) to the full set of B° D' K°s-k-)-k+ candi- 
dates. The dashed line is the fitted background PDF while the 
solid continuous line is the sum of the signal and background 
PDFs. The solid dots are data used in the fit. 



To extract the branching fractions for the decay modes 
^0 ^ K*+TT- and D-{-^ K°sTt-)tt+ the following 

equation was used: 




2.0 3.0 
Fisher 



1.0 2.0 3.0 
Fisher 



B = 



Nbb X e ' 



(3) 



FIG. 2: Maximum likelihood fit projections of niES, Ai5 and 
for signal enhanced samples of charmless B° — > K'^7v'^-k~ 
and B° K*^tv~ candidates. The dashed line is the fitted 
background PDF while the solid continuous line is the sum 
of the signal and background PDFs. The solid dots are data. 



where n^ig is the number of signal events fitted, e is the The left column (plots a-c) has the B° if^Tr+yr" projec- 
tions and the right column (plots d-f) has the — » K*^tt~ 
projections, the top, middle and bottom rows being the tues, 
AE and J- distributions respectively. 



signal efficiency obtained from MC and Nbb is the to- 
tal number of BB events. For the charmless inclusive 
B'^ — > K^TT^TT~ branching fraction, the efficiency varies 
over the Dalitz plane and the distribution of events across 
it is a priori unknown, consequently the total efficiency 
is unknown. Therefore, to calculate the branching frac- 
tion, a weight was assigned to each event such that, for 
the jth event Wj = Y^^Vsig,iPi{d,Xj)/Y,^nkPk{a,Xj) 
where Vsig^i is the signal row of the covariance matrix 
obtained from the fit This procedure is effectively 

a background subtraction where these weights have the 
property J2j — ''^sig ■ The branching fraction is then 
calculated as B = J2j'^j/i^j ^ ^bb) where ej is the 
efficiency which varies across the Dalitz plot and is sim- 
ulated in small bins using high statistics MC. 

Figure |21 shows the fitted projections for both 
K^TT+TT" and B° K*+n~ candidates, whilst the fitted 
signal yield and measured branching fraction are shown 
in Table HI for all the modes under study. The average 
efficiency for —f K'^tt^tt" signal events was approxi- 
mately 8%. Figure O shows the signal mass projections 
of rrii^o^ using — > K'^tt^tt^ candidates. The m^^o^ 
distribution clearly shows a peak at 0.9 GeV/c^ which 
corresponds to the K*^{892) and there is a broad struc- 



TABLE I; Signal yields and branching fractions for — > 
A"°7r+7r-, B° K'+n- and B° D"(^ K'^tv-)-jt+ where 
the first error is statistical and in the case of the branch- 
ing fraction measurements the second error is systematic. 
The efficiency of selecting B° K*+n- and B° D' {-^ 
K%-R~)-R^ events was found to be 5.1% and 12.4% respec- 
tively. The B^ — » K*^'K~ branching fraction takes into ac- 



count that B{K*+ 
metry. 



K^-K^) — 2/3, assuming isospin sym- 



Mode 


Signal Events 
Yield 


Branching Fraction 
(x 10-'') 


B° ^ K'+^- 
B° D-{-* KIti-)ti+ 


310 ± 27 
59 ± 11 
472 ± 24 


43.7 ± 3.8 ± 3.4 
12.9 ± 2.4 ± 1.4 
42.7 ± 2.1 ± 2.2 



Contributions to the branching fraction systematic er- 
ror are shown in Table ^ Errors due to pion tracking, 
particle identification and K^^ reconstruction efficiency 
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FIG. 3: The rriKTr distribution of _B" K^\^tv~ candidates, 
weighted by W such that background events are subtracted. 
The one-dimensional distribution is obtained by considering 
events on a Dalitz plane with axes rn^j^o^^ and ?n^o^_. The 

^2 



two axes are merged into one (mjf^) by folding the Dahtz 
plane along the line corresponding to m^o^^ ~ "i-xOtt- 
order to obtain the above niKn mass distribution. 



were assigned by comparing control channels in MC and 
data. To calculate errors due to the fit procedure, a large 
number of MC samples containing the amounts of signal, 
continuum and i?-background events measured or fixed 
in data were used. The differences between the generated 
and fitted values using these samples were used to ascer- 
tain the sizes of any biases. Small biases of the order of 
a few percent were observed that were a consequence of 
small correlations between fit variables and were there- 
fore assigned as systematic errors. The uncertainty of 
the B-background contribution to the fit was estimated 
by varying the measured branching fractions within their 
errors. Each background was varied individually and the 
effect on the fitted signal yield was added as a contribu- 
tion to the uncertainty. For i?° K*^Tr~ there was also 
the S-background contributions from higher kaon reso- 
nances which was obtained from fits to data and added 
as a systematic. The uncertainty due to simulated PDFs 
was obtained from the channel 5° D^{-^ K^tt^)-k^ 
and by varying the PDFs according to the precision of 
the parameters obtained from MC. In order to take cor- 
relations between parameters into account, the full cor- 
relation matrix was used when varying parameters. All 
PDF parameters that were originally fixed in the fit were 
then varied in turn and each difference from the nominal 
fit was combined and taken as a systematic contribution. 
The error in the efficiency was due to limited MC statis- 
tics, where over a million MC events were generated for 
the decay B° K° TT^TT and over one hundred and fifty 
thousand MC events were generated for the decays 
K*+TT- and 5° D-{-^ i4:°7r-)7r+. The same uncer- 
tainty due to the error in the number of BB events was 
added to all channels. 



TABLE IL Summary of systematic uncertainty contributions 
to the branching fraction measurements — > K^-k'^'k~ , i?" 
— > K*'^-K~ and 5° — > > K^ti~)-k'^ . The errors are shown 

as a percentage of the measured branching fraction. 



Error 


— » KO-k+-it~ 


B" — > K*+-K~ 


5° — > D~TT + 


source 


error (%) 


error (%) 


error (%) 


Tracking 


1.7 


1.7 


1.7 


Particle ID 


1.9 


1.2 


3.0 


Efficiency 


4.2 


3.5 


3.0 


Fit Bias 


4.1 


3.3 


1.8 


B-background 


3.6 


9.0 


0.3 


PDF params. 


1.5 


0.5 


0.4 


Efficiency 


1.7 


0.6 


0.6 


No. of BB 


1.1 


1.1 


1.1 


Total 


7.8 


10.5 


5.1 



Interference was also considered for the decay B^ 
K*+Tr~ where effects between the K*'^{892) and S wave 
final states (non-resonant and Kq^{14:30)) cancel and 
the K*+{892) and D wave final states (iC2 + (1430)) can- 
cel. This is not the case for P wave amplitudes such as 
A'*+(1410), yet this effect was considered negligible due 
to the small branching fraction of _R'^~''(1410) — > K'^n'^ 
(6.6 ± 1.3% 13). 

The CP violating charge asymmetry for the decay 5° 
^ K*+TT~ was measured to be Ak-t, = 0.23 ± 0.18lJ];[!^, 
where the first error is statistical and the second errors 
are systematic. The background asymmetry A^'l^^ was 
measured to be 0.01 ± 0.01 and as a further study the 
asymmetry Adtt for D~(— » K'^Tr~)7r~^ was mea- 
sured to be 0.00 ± 0.05 and the background asymmetry 
■^Lhr' 0-^^ ^ 0.04, were the errors are statistical only. 

The systematic error on Ak*tt was calculated by con- 
sidering contributions due to track finding, particle iden- 
tification, fit biases and i3-background asymmetry uncer- 
tainties. Biases due to track finding and particle identi- 
fication were found to be negligible. The fit bias con- 
tribution to the systematic error was calculated using 
a large number of MC samples. The contribution from 
-B-background was calculated by varying the number of 
expected events within errors and by assuming a conser- 
vative CP violating asymmetry of ± 0.5 as there are no 
available measurements for these decays. The resulting 
systematic uncertainty on the asymmetry was measured 

to be +° °^ 
lo ue _Q 06- 

In summary, the branching fractions for B^ — )■ 
K^TT+TT-, B° K*+Tr- and 5° D-{~^ K'^^tt-)tt+ 
decaying to a K'^tt^tt^ state have been measured and 
agree with previous measurements 0, S • The direct 
CP violating parameter Ak^it was measured for the de- 
cay — > K*'^w~ and is in agreement with the CLEO 
measurement 0, with no evidence of CP violation with 
the statistics used. Using larger datasets, one can extract 
amplitudes and relative phases of the resonant contribu- 
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tions to the Dalitz plot, with the possibiHty to observe 
new B-meson decays. 
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